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ABSTRACT 


Two  methods  for  the  removal  in  real  time  of  stable  spectral 
lines  from  seismic  data  were  investigated.  The  first  method  involved  the  gen¬ 
eration  of  a  cosine  wave,  by  means  of  a  digital  feedback  sy^em,  to  approximate 
a  spectral  line.  In  this  method,  the  cosine  genera  or  is  s  ijusted  in  accordance 
with  a  mean- square -error  criterion.  This  system  of  adjustnoviit  was  simulated 
and  proved  to  be  unstable.  In  the  second  method,  a  Widrow  adaptive  prediction 
filter  is  used  to  remove  the  deterministic  component.  A  simulation  of  this  pre¬ 
diction  method  was  carried  out  using  one  seismic  data  channel  from  a  short- 
period  array.  Results  show  that  this  method  significantly  attenuates  some  of 
the  spectral  lines.  However,  if  a  signal  were  present,  some  amplitude  and 
phase  distortion  would  be  caused  in  the  signal.  In  addition,  this  system  would 
be  cumbersome  to  implement. 
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SECTION  I 

INTRODUCTION  AND  SUMMARY 


This  report  discusses  the  removal  of  stable  spectral  lines  from 
seismic  data  in  real  time.  The  source  of  these  lines  is  suspected  to  be  mining, 
drilling,  or  other  man-made  activity.  At  some  stations,  these  spectral  lines 
are  of  relatively  high  energy  and  are  present  in  the  band  of  strong  signal  energy. 
Because  the  exact  frequency  of  each  line  is  not  known,  and  because  these  fre¬ 
quencies  may  slowly  fluctuate,  a  single  fixed  filtei  with  appropriately  positioned 
nulls  may  not  be  effective.  Furthermore,  it  would  be  difficult  to  implement 
such  a  fixed  filter  to  obtain  a  relatively  distortion-free  response  over  the 
remaining  spectrum. 

A  time-domain  technique  for  predicting  extremely  narrow  spec¬ 
tral  line  components,  so  that  they  mav  be  removed  by  subtraction,  has  been 

if 

demonstrated  in  a  previous  report.  This  technique  has  been  shown  to  be  highly 
effective,  and  its  use  would  be  most  beneficial  at  stations  having  narrow  spec¬ 
tral  lines  in  the  ambient  noise  at  frequencies  in  the  signal  band.  However,  this 
technique  involved  designing  an  optimum  100-point  linear  prediction  filter  from 
over  2000  points  of  data  and  would  be  relatively  difficult  to  implement.  The 
problem  investigated  in  this  report  is  to  find  a  simplified  method  which  would 
be  suitable  for  real-time  processing.  If  such  a  method  works,  the  spectral 
lines  in  the  output  data  should  be  greatly  attenuated  in  comparison  to  the  input 
data,  and  the  rest  of  the  frequency  spectrum  should  be  relatively  undistorted. 

Two  methods  for  the  removal  of  spectral  lines  were  studied. 
Assumed  for  each  method  are 


if 
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•  The  data  arc  the  sum  of  two  processes:  a  random, 
wide -band  process  and  a  relatively  deterministic 
process 

•  The  relatively  deterministic  process  is  that  which 
gives  rise  to  the  spectral  lines  and  can  thus  be 
predicted  and  removed  from  the  data 

A  digital  feedback  system,  called  a  cosine  generator,  is  used 
in  the  first  method  to  generate  a  cosine  wave.  The  purpose  is  to  approximate 
the  spectral  line  by  a  generated  cosine  wave,  and  the  cosine  generator  is  ad¬ 
justed  in  accordance  with  a  mean- square -error  criterion.  This  system  of 
adjustment  has  been  simulated  and  proves  to  be  unstable. 

In  the  second  method,  a  Widrow  adaptive  prediction  filter  is  used 
to  remove  the  deterministic  component.  A  simulation  of  this  prediction  method 
has  been  carried  out  using  one  seismic  data  channel  from  a  short-period  array. 
The  results  show  that  this  method  does  attenuate  significantly  some  of  the  spec¬ 
tral  lines.  However,  some  amplitude  and  phase  distortion  would  be  caused  in 
a  signal,  if  it  were  present.  In  addition,  it  would  be  cumbersome  to  imple¬ 
ment  this  system. 


SECTION  II 

METHOD  1:  COSINE  GENERATOR,  MEAN-SQUARE -ERROR  CRITERION 
A.  GENERAL  DESCRIPTION 

In  this  method,  a  digital  feedback  system  is  used  to  generate 
a  sinusoid  which  is  used  to  approximate  a  spectral  line  in  the  data.  The 
error  of  the  approximation  is  calculated  and  statistical  parameters  of  the 
signal  are  estimated  at  regular  intervals.  These  quantities  are  employed 
using  a  heuristic  least  mean-square-error  criterion  for  adjusting  the  cosine 
generator.  The  scheme  is  heuristic  in  that  certain  functional  dependences 
are  neglected  for  simplicity. 


Actual  simulation  shows  this  type  of  system  to  be  unstable.  It 
is  expected  that  the  neglected  functional  dependences  cause  this  instability. 
Although  these  dependences  were  obvious,  it  was  hoped  that  their  presence 
would  not  adversely  affect  the  proposed  system. 

B.  ANALYSIS  OF  OPERATION 


Let  x^  be  the  input  data,  and  r  be  the  generated  sinusoid. 


Then  rk  is  given  by  r^  =  B  r^  -  A  ru  ,.  The  error  in 


k-2 


approximating  x 


With  rk  is  \  “  *k  “  rk*  The  squared  error  is  given  by 
ek  -  [>•  -B>  A] 
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If  the  sequence  were  not  f  u>c£ionally  dependent  on  A,  B,  the 
quantities  in  the  matrix  of  Equation  2-2  would  not  be  determined  by  A,E.  If 
this  is  assumed,  A  and  B  can  be  considered  undetermined  constants,  and  the 
mean- squared-error  can  be  minimizesd  by  selecting 


-B 

A 


2  1 

r  r  r. 

x  r 

k-1  k-1  k-2 

k  k-1 

r.  .  r,  „  r^ 

x  r 

k-1  k-2  k-2 

k  k-2 

(2-3) 


Estimates  of  the  covariances  must  be  made  to  implement  this  system,  and  thus 
the  system  error  and  stability  are  sensitive  to  how  these  estimates  are  made. 
The  covariance  matrix  of  Equation  2-2  is  named  R,  and  an  estimate  of  it  is 
made  at  regular  intervals  in  accordance  with  the  following  relationship. 
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0<e<  1 


This  method  of  estimation  is  effectively  the  same  as  computing  the  covariance 
matrix  using  data  which  is  weighted  by  a  decaying  exponential  function.  The 

A 

initial  estimate  R  is  selected  to  be  the  identity  matrix. 


The  oscillator  is  initialized  to  generate  a  sine  wave.  This  is 
done  by  selecting  rQ  =  0,  r^  =  1,  and  A  =  1.  The  coefficient  B  is  set  by  the 
relationship 


1 

2B 


COS  Cp 


,  where  Cp  = 


tt  frequency  °* 

folding  frequency 


In  this  relationship,  the  folding  frequency  is  considered  to  be  one-half  the 
sampling  frequency. 
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SECTION  in 

METHOD  2s  WIDROW  PREDICTION  FILTER 

This  method  uses  a  Widrow  prediction  filter  to  predict  the  sum 
of  aU  the  spectral  lines  in  the  data  *  The  filter  is  applied  to  the  data,  and  the 
system  output  is  formed  by  subtracting  the  prediction  of  the  spectral  lines  from 
the  input  data.  This  output  is  thus  constructed  to  approximate  the  input  data 
minus  the  spectral  lines.  Because  this  output  is  the  "error"  between  the  pre¬ 
diction  and  actual  data,  this  system  is  sometimes  referred  to  as  a  "prediction- 
error  filter". 


A  simulation  of  this  prediction  method  for  spectral  line  removal 
has  been  carried  out  using  seismic  data  from  a  short-period  array.  The  input 
data  spectrum  is  shown  in  Figure  HI-l.  A  20-point  long  prediction  filter  was 
adaptively  designed  to  predict  the  data  several  points  in  the  future.  In  applying 
the  filter,  several  experiments  were  done  using  various  prediction  lengths. 

(The  design  parameter  determining  which  point  in  the  future  to  predict  is  re¬ 
ferred  to  as  the  prediction  length.  For  example,  the  20-point  filter  operating 
on  data  points  1  through  20,  and  designed  to  predict  data  point  25,  would  have 
a  prediction  length  of  five.  In  operation,  the  error  in  prediction  would  be  used 
to  update  the  filter,  and  the  filter  would  be  applied  to  points  2  through  21*  and 
the  output  would  be  the  prediction  of  point  26. )  Also,  the  rate  of  adaption, 
referred  to  as  the  convergence  parameter,  was  changed,  yielding  a  table  of 
results.  Specifically,  the  prediction  lengths  were  10,  50,  100,  and  200  sample- 

points  ahead  and  the  filter  convergence  parameter  k  took  cn  values  of  0.  005 

s 

Xkmax*  ^^Xkmax'  ^^X^max'  where  ^max  is  a  theoretical  limitation 
on  the  convergence  factor.  The  output  spectra  are  shown  in  Figures  HI-2 
through  III— 1 2. 

It  is  seen  from  Figure  III — 1  that  by  far  the  greatest  proportion 
of  energy  from  spectral  lines  occurs  in  those  two  lines  between  2.  5  and  3.  0  He. 

£  " 

Widrow,  Bernard,  1966:  Adaptive  Filter  s  I:  Fundamental,  Stanford  University 
Tech.  Rpt.  No.  6764-6,  Systems  Theory  Laboratory,  Contract  No.  NOBsr-95038 
and  Contract  No.  DA-01 -021  AMC-9001 5(Y),  Dec. 


Figure  III-6  shows  that  these  lines  are  attenuated  by  13  db  when  the  data  are 

filtered  with  a  prediction-error  filter  which  predicts  50  points  ahead  with  a 

convergence  parameter  equal  to  0.  05 xk  .  While  other  spectral  lines  are 

max 

not  significantly  attenuated ,  it  is  noted  that  these  lines  contain  relatively  little 
energy  in  comparison  to  the  larger  two  lines  anti  the  results  were  obtained  on 
a  comparatively  short  piece  of  data. 

The  data  are  3300  points  long  (about  4  min)  a.'vd  the  method  used 

to  filter  it  follows.  The  initial  filter  is  chosen  to  be  composed  only  of  zeros. 

As  the  filter  processes  the  first  1000  points  of  data,  the  convergence  parameter 

kg  decreased  linearly  to  its  final  value.  This  new  value  is  kept  constant  for 

the  remainder  of  the  processing  in  that  trial.  (The  initial  values  of  k  depend 

upon  the  corresponding  final  values.  For  final  values  of  0.  05  k  and  0.  25 

max 

^max'  value  of  k^  is  twice  these  corresponding  values.  For  the 

final  value,  0.  005  k  ,  k  was  initially  set  to  0.  5  k  . ) 

max  s  7  max 

To  determine  how  much  the  prediction  filter  would  distort  a  sig¬ 
nal  if  a  signal  were  present,  the  frequency  response  of  several  prediction  filters 
were  computed.  Because  the  filters  are  adaptive,  they  change  as  they  piocsss 
each  new  point  of  data,  and  the  terminal  filters  were  selected  for  computing 
the  frequency  response.  The  responses  appeared  similar,  and  the  response  of 
the  Biter  which  produced  the  output  spectrum  shown  in  Figure  HI -6  is  shown  in 
Figure  III- 13.  The  relatively  large  amplitude  response  throughout  the  frequency 
spectrum,  combined  with  a  nonuniform  phase  response  (not  shown),  indicate  a 
considerable  amount  of  distortion  would  occur  in  a  signal  if  it  were  present. 

In  summary,  the  results  do  not  indicate  that  all  of  the  spectral 
lines  are  significantly  attenuated.  In  addition,  amplitude  and  phase  distortion 
would  be  caused  in  a  signal  if  it  were  present.  Although  results  from  processing 
a  longer  segment  of  data  may  show  a  more  uniform  attenuation  of  all  spectral 
line 8,  the  problem  of  signal  distortion  would  still  be  present. 
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Figure  III- 5.  Output  Power  Spoctrum,  kg  =  0.  005,  Prediction  Length  =  50 
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Figure  III-6.  Output  Power  Spectrum,  kg  =  0.  05,  Prediction  Length  =  50 
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SECTION  IV 

CONCLUSIONS  AND  RECOMMENDATIONS 

Neither  of  the  methods  examined  has  proved  to  be  satisfactory  in 
the  extraction  of  spectral  lines  from  seismic  data.  A  simple  method  for  real¬ 
time  elimination  of  spectral  lines  would  probably  require  a  significant  investi¬ 
gation.  Since  that  is  beyond  the  scope  of  the  very  small  effort  assigned  to  this 
task  tinder  the  present  contract,  it  is  recommended  that  the  study  of  spectral 
line  extraction  be  discontinued. 
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eration  of  a  cosine  wave,  by  means  of  a  digital  feedback  system,  to  approx¬ 
imate  the  spectral  lines.  In  this  method,  the  cosine  generate"  is  adjusted 
in  accordance  with  a  mean- square -error  criterion.  This  system  of  adjust¬ 
ment  was  simulated  and  proved  to  be  unstable.  In  the  second  method,  a 
Widrow  adaptive  prediction  filter  is  used  to  remove  the  deterministic  com¬ 
ponent.  A  simulation  of  this  prediction  method  was  carried  out  using  one 
seismic  data  channel  from  a  short-period  array.!  He  suits  show  that  this 
method  significantly  attenuated  some  of  the  spectral  lines.  However,  if 
a  signal  were  present,  some  amplitude  and  phase  distortion  would  be  caused 
in  the  signal.  In  addition,  thi*  system  would  be  cumbersome  to  implement. 
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